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Abstract 



In this paper we present an exact study of the relaxation dynamics of the backgammon 
model. This is a model of a gas of particles in a discrete space which presents glassy 



^ ■ phenomena as a result of entropy barriers in configuration space. 



The model is simple enough to allow for a complete analytical treatment of the dynamics 
in infinite dimensions. We first derive a closed equation describing the evolution of the occu- 
pation number probabilities, then we generalize the analysis to the study the autocorrelation 
function. 

We also consider possible variants of the model which allow to study the effect of energy 
barriers. 
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I. INTRODUCTION 



The nature of the glass transition is still poorly understood []T],@]- Under slow cooling 
real glasses reach a metastable phase of free energy larger than that of the crystal phase. 
Glasses show a strong slowing down of the dynamics when the temperature is lowered 
and the transport coefficients increase by several orders of magnitude in a narrow range 
of temperatures. It is natural to think that the appareance of high free-energy barriers is 
the mechanism responsible for the glass transition. But free-energy barriers are composed 
of energy barriers and entropy barriers. The question about the relevance of both kind of 
barriers in real glasses is of the outmost importance. Activated jumps of energy-barriers 
are strongly dependent on temperature. The typical time r to overcome an energy barrier 
AE is r ~ exp(^P) where T is the temperature. This typical time diverges when the 
temperature T goes to zero. Conversely, relaxation times related to entropy barriers do not 
depend directly on the temperature. 

The simplest way to visualize entropy barriers is the following, Consider a dynamics in 
which at each time step the system can reach a new state with uniform probability; the 
typical time to decrease the energy of one unity is r ~ ^f- = expfAS) where AS* is the 
height of the entropy barrier stands for the initial available volume of phase space and 
Qf stands for the final volume of phase space with lower energy). 

If the effect of energy and entropy barriers is combined one expects that entropy barriers 
should affect the temperature activated relaxation time in its prefactor r ~ ^-exp(^). 
According to that, the relaxation time can diverge independently of the temperature if the 
phase space volume of lower energy configurations in the system shrinks to zero during the 
dynamical evolution. The idea that an entropy crisis could be relevant to the glassy transi- 
tion is very old HU, and it has had interesting developements in recent times |5]|| in the 
framework of mean-field theory of disordered systems. However in the models studied in 
HH it is very difficult to disentangle entropic effects from energetic ones. To this aim a 
simple model (the backgammon -BG- model) was recently proposed by one of us || (here- 
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after referred as I), in which energy barriers are completely absent (a diffusive model with 
entropy barriers has also been considered in |ID|). While the model has no thermodynamic 
transition, it shows a slow dynamics with strong hysteresis effects and Arrhenius behavior of 
the relaxation time. The off-equilibrium dynamics of this model was studied subsequently 
by us |J (hereafter referred as II) using an adiabatic approximation, obtaining fair good 
results concerning the relaxation of the energy. The same approximation has been recently 
rederived, and slightly refined, in a paper by Bouchaud, Godreche and Mezard |JTTJ . 



In this paper we derive the exact mean-field equation for the order parameter for the 
dynamics of the BG model, which turns out to be the energy itself. The techniques we use 



are similar to these of ||11|| , however, the equations we get were not discussed there. We find 
that the energy verifies a causal functional equation with memory. This is at variance with 
the approximate treatments where the evolution is described by a Markovian equation. 

In its original formulation, the model does not have any energy barriers. However, in 
real systems energy barriers are present. The BG model is flexible enough to allow for the 
introduction (and the tuning) of energy barriers. This is done by simple modifications of 
the Hamiltonian of the system. The formalism developed for the original model applies in 
these cases. 

In the second section we define the Hamiltonian of the BG model and the Monte Carlo 
dynamics we have used to study it. In the third section we present some exact results for the 
behavior of the one-time quantities (for instance, the energy) and for the two-time quantities 
like the correlation function. The fourth section is devoted to the study of the effect of energy 
barriers in the BG model. Finally we present the conclusions and a discussion of the results. 

II. THE BG MODEL AND THE DYNAMICS 

Let us take iV distinguishable particles which can occupy M different states and let us 
denote by p = the density, i.e. the number of particles per state. The BG Hamiltonian 
is defined by, 
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M 

H=~Y. S nr,0 (1) 
r=l 

where n r is the occupation level of the state r = 1, ...,M, i.e. the number of particles 
which occupy that state. The numbers n r satisfy the global constraint, 

M 

J2n r = N . (2) 

r=l 

Eq.([l|) shows that energy is simply given by the number of empty states (with negative 
sign). We define the occupation probabilities, 

1 M 

p k = m E < K,k > (3) 

which is the probability of finding one state occupied by k particles. The statics of this 
model in the canonical ensemble can be easily solved (see (I) and (II)). In particular one 
gets the result, 

p k = z*- 1 exp(p6 kfi ) 
' k\ exp(z) 

where z is the fugacity and j3 is the inverse of the temperature T and they are related 
by the condition, 

p{e P - 1) = (z - p)e z . (5) 

expressing that the density is fixed to p. 

The probabilities Pk satisfy the relation Pk — 1 an d they yield all the static ob- 
servables, in particular the energy U = —Pq. Several dynamical rules, thermalizing to 
the Boltzmann distribution, can be attached to the model. The simplest choice (I) is the 
Metropolis single particle dynamics, in which at each sweep a particle is chosen at random 
and it is proposed a move to a new state. The move is accepted with probability one if the 
energy does not increase and with probability exp(— (3) otherwise. 

In the mean-field version of the model, the possible arrival states of the particles are 
chosen at random with uniform probability in all the space. This random motion of the 



particles allows a complete analytical treatment of the problem []. Finite dimensional models, 
where at each sweep the particles are only allowed to move to neighbours on a lattice are 
currently under study J7|. 

The model has no energy barriers. Consequently there is no frustration (in the usual 
sense) and no metastable states. However it was shown in (I) that the dynamics is highly 
non trivial and a dramatic slowing down occurs at low temperatures. This can be qualita- 
tively understood as follows. Suppose the system is at zero temperature and the dynamics 
starts from a random initial configuration of high energy. As the system evolves towards 
the equilibrium more and more states are progressively emptied and the energy decreases. 
Because the average number of particles per occupied state increases with time (the total 
number of particles is conserved) the time needed to empty one more state also increases. 
The result is that the energy goes extremely slowly to its equilibrium value. 

The dynamical quantities we are interested in are the time-dependent occupation number 
probabilities 

i M 

= 77 E <**(«),*> (6) 
m r =l 

(E(t) = —Po(t)) and the two time energy-energy correlation function ||, 

C E (t, s) 

t>s (7) 



Er $n r (t),o8n r (s),0 ~ E(t)E(s) 

-E(s)(l + E(s)) 
P{n r {t) = 0,n r {s) = 0) - P (t)P (s) 



P (s)[l - P (s)} 

At finite temperature, when t,s » t eg ~ exp(/3)//3 2 (see reference (II) and also section 
III.B) this function is time translationally invariant. In the regime in which both times t, s 
are much less than t eq , and at all times at zero temperature, the system is off-equilibrium, 
time-translation invariance does not hold, and the correlation function displays aging (see 



The interesting case of a sequential dynamics is more complicated. 
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III. MEAN-FIELD EQUATIONS FOR THE DYNAMICS OF THE BG MODEL 



In this section we derive exact mean-field equations for the Monte Carlo dynamics of the 
BG model. First we adress the dynamical problem associated to the one-time probability 
distributions Pk{t). These probabilities generate an infinite hierarchy of Markovian equations 
which can be closed in terms of the only quantity Po(t). Then we will study the two- 
time correlation functions in a similar way. For simplicity, we will restrict all the future 
computations to the case p = 1 (i.e M = N), the generalization to an arbitrary density 
being very simple. 



A. Dynamical equations for P k (t) 

The purpose of this section is to write the dynamical evolution equations for the prob- 
abilities Pfc(t) and, in particular, for the internal energy E(t) = — Pq(£). An elementary 
Monte Carlo move consists in a random selection of one particle (hence, the probability to 
select a particular departure state d is n^/N where is the occupation level of that state) 
and moving it to a randomly selected arrival state a with uniform probability independent of 
the occupation level n a . One Monte carlo step (our unity of time) consists of N elementary 
moves. In the elementary move there are several processes which contribute to the variation 
of Pfc(t). In appendix A we write explicitly the balance equations, the result is: 

dP (t) 

— ^ = (k + l)(P fc+1 - P k ) + P fc _x + P (e^ - 1)(5 M - 4,o - kP k + (k + l)P fe+1 ) (8) 

where the time index for the probabilities P k has been omitted. The previous equation holds 
for k > with P_i = 0. In particular for k = we obtain the equation studied in (II), 

dP 



Of 



P 1 (l-Po)-e^P (l-Pi) (9) 



The hierarchy was closed in (II) assuming fast relaxation on the surfaces of constant energy, 
and slow variation of the energy itself. In this condition eq.@ was solved assuming P k {t) = 
exp[/9(t)4,o — z W\^h. — with f3{t) and z(t) related at all times by eq.@. 
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Here we study the hyerarchy @ with the method of the generating function, that we 
define as 

oo 

G(x,t) = Y,x k P k (t) (10) 

fc=0 

A similar approach was also used in [|TTJ where the adiabatic approximation of (II) was 
rederived and improved?]. 

From the equation (|8|) it is easy to check that the G(x, t) satisfies the partial differential 
equation, 

= {x _ 1)[G(a . t) + A(t) _ (1 + \{t))™jg*>] (11) 

with X(t) = Po(t)(e~P — 1). Eq. (|Tl"l) is a non linear partial differential equation, the 
nonlinearity is contained in the dependence of A on Po(t) = G(0,t). 

The equilibrium solution G eg (x) is easily obtained from equations @ and (|T0l), 



P P _ i _i_ P zx 

G eq (x) = J/ (12) 



and one can check that this is consistent with eq. ([TT|) . 

The previous partial differential equation can be implicitly solved to get G(x,t) as a 
functional of A. The details are presented in the Appendix B, we give here the result 

G(x, t) = e (a; - 1)D(t ' 0) G (l + (x- l)B(t, 0)) + (x - 1) f ds\(s)B(t } s) e ( "- 1)D(t ' s) (13) 

Jo 

where we have written 



B{t, s) = exp (- J* dv [1 + X(v)} 
D(t,s) = f dv B(t,v) (14) 

J s 



and G (x) = G(x, 0) is the initial condition at time t — 0. Setting x = in ( p~3|) we get 
a closed equation for Po(t) which reads 



2 The technique of the generating function in the study of the dynamics has also been applied to 
some mean- field spin glass models [12]. 
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P (t) = e- D ^G Q (l - B(t, 0)) + (1 - e'P) f dsP (s)B(t, s) e D ^ (15) 

J 

The previous equation, although strongly non-markovian is causal, as the l.h.s. depends 
on the values of Pq(s) for s < t. It has a unique solution that can be found numerically 
with good precision, discretizing the time and integrating it step by step. The evaluation 
of the previous expressions gives the full solution of the BG model as far as the one-time 
dynamical quantities are concerned. 

The solution of (|13|) is explicit at infinite temperature (/3 = 0). In this case \(t) = and 
the solution of G(x,t) simplifies, 

G(x, t) = eCi-OO-i) G ((x - l)e-* + 1) (16) 

It is not surprising that at infinite temperature the system goes exponentially fast to the 
equilibrium (with relaxation time equal to 1). At infinite temperature the equilibrium prob- 
abilities eq.([|) are given by Pk = tj-, the energy being E = —P — — K If we start from the 
initial condition in which all particles occupy the same state (Pq — 1, = 0, k > 0) then 
we have Gq(x) = 1. From eq.fllBD we obtain the time evolution of the energy, 

E(t) = -G{0,t) = -e 6 "- 1 (17) 

We studied numerically the solution of (|1|) at T = 0. In figure 1 we display the result for 
the energy, starting from the initial condition Pfc(0) = l/(ek\) at time (i.e. Gq(x) = e x ~ l ). 
For comparison we plot the results of the Monte Carlo simulations and of the adiabatic 
hypothesis of (II) with the same initial condition.^ 

B. The Correlation Function Cg(i, s) 

In this section we investigate the behavior of the energy-energy correlation functions (^). 
We proceed in a similar way as we have done for the occupation probabilities. We need to 



3 The adiabatic hypothesis gives better results if the integration is started at latter times. 
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study the joint occupation probability in a given site r at two different times t, s (t > s), 
P(n r (t) = 0,n r (s) = 0) = P{n r {t) = 0|n r (s) = 0)P (s). The correlation function eq. (|7|) 
can be written as 

r , P(n T .(t) = OMg) = 0) - Po{t) 

C E (t, s) = — ^ . (18) 

We now write a set of equations that allow to study P(n r (t) = 0\n r (s) = 0). 
Let us define the probabilities 

u k (t, s) = P(n r (t) = k\n r (s) = 0) (19) 

i.e. the occupation number probabilities in the set S s of states which are empty at time s. 
In general, it is possible to restrict the balance equations that led to fl8|) to any subset S of 
the whole space. Irrespectively of S the result is: 

^ = - u k + [{k + l)u k+1 - kv k \[l - P (l - e- 13 )] 

-(4,i - 4,o)h(l - Pi) + ViP ](l ~ e-* 3 )- (20) 

In particular if the set S is the whole space v k = P k and we get back to 

Of course the initial conditions depend on the set under study. For the set S s we are 
interested to, we must choose 

u k (s, s) = 8 k}0 . (21) 
In terms of the generating function 

oo 

r(x,M) = ]>> fc z/ fc (M) (22) 

fc=0 

eq. (^0|) reads 

Br 8r 

— = (x-l)[T-(l- P (l - e-?))— - (Ml ~ Pi) + ^o)(l - e-?)] (23) 
with condition at time s 

F s (x) =T{x,s,s) = 1. (24) 



Note that if we suppose to know the Pk(t) then the eg. ( pO| , p3D are linear. Obviously, if one 
considers the set S s complementary to S s , and its respective generating function T, it holds 
the equality: P (s)T(x, t, s) + (1 — P (s))T(x, t, s) = G(x, t, s). 
Defining 

H{t, s) = h(t, s)(l - Pi(t)) + vt{t, s)P (t)](l - e~ p ) 
Bit, s) = exp (- J* dv (1 - P (*)(l - e^))) 
D(t,s) = f dv B(t,v) 



(25) 



we find 



F(x, t, s) = e {x ~ 1)D{t ' s) - (x - 1) / dufj,(u, s)B{t, u) e^ D ^ 



(26) 

which depends implicitly on v and V\. In order to find a closed system we have to consider 
eq. (|26"D and its x-derivative in x = 



v (t,s) 
vi(t, s) 



1+ C dv 

t 



-P (U, S)[l - (1 - Pl{u)){\ - + S) 



du 



fi{u, s)B(t, v)e- D{t ' u) {D{t, u) - 1) + e~ D ^D(t, s 



(27) 



The system (p7D , if one assumes known the probabilities Pk{t), consists in a vectorial 
linear Volterra equation of second kind for vq and V\ which can in all generality be integrated 
numerically, and in some particular case also analytically. 

The simplest case is the equilibrium at finite temperature. In that case, = P^ 9 and 
the various functions appearing in (137]) are time traslation invariant. In these conditions 
eq.(^TJ) can be solved in Laplace transform. Simple algebra, and the formula (see e.g. [O]) 



/ dt exp(— aexp(i) — E) = a~ E j(p, a) 
Jo 



(28) 



(7 is the incomplete gamma function), shows that uq(E), the Laplace transform of u (t — s 
is given by 

A(E) + 2-1(1 - EA{E)) 



ME) 



(z-l)e z 



(z-l)e z +l 



+ ^E (l-EA(E)) 



(29) 
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where we have expressed all the equilibrium quantities in terms of the fugacity z (see eq. 
(D), E is the Laplace variable conjugated to time, and 

A(E) = — \—(\^e\ (30) 



e z z zE-l j Q 



v (E), as it should, has a pole in E = with residue P coresponding to vo(t) — > P f° r large 
time. Poles on the real negative E axis correspond to exponential relaxation modes. The 
largest relaxation time is given by minus the inverse of the value of E in the pole closest to 
the origin. This can be obtained explicitely for large /3, where z ~ (3 — log(/3) is large, from 
the asympotic expantion of EA(E) for small E 

EA(E) « e- z + E. (31) 

The result is simply E po i e ~ — e~ z and correspondingly r max rs e z ~ exp{(5)/ (5. 

In the off-equilibrium regime the integration of ( |2T| ) can be performed numerically. In fig. 
2 we show the result of the integration for T = for different values of s ( i.e. different waiting 
times) compared to the Monte Carlo results. Although we did not try very sophisticated 
algorithms, with standard ones we were able to reach enough large times to see the 



scaling behavior Cs(t, s) = f((t — s)/s) observed numerically in I. It would be interesting 
to see if equation (Pq ) could be solved with the aid of some simple approximation as it is 
the case for the energy (II and 01 ) • 



IV. THE EFFECT OF ENERGY BARRIERS 

The BG model has no energy barriers and hence there is no finite temperature thermo- 
dynamic phase transition. In real glasses energy barriers are always present and it can be 
instructive to understand their effect when combined with entropy barriers. One can easily 
modify the Hamiltonian ([]]) to include energy barriers. In this paper we have focused on 
two different ways. In the first, we have considered interaction between the different states, 
introducing an energy gain when groups of states are simultaneously empty. This interac- 
tion term is enough to make appear a finite temperature thermodynamic transition, but 
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metastability and frustration are absent and the system monotonically reaches the ground 
state at zero temperature. In the other, we introduced metastable configurations in the 
dynamics. In this case the system fails to reach the ground state at zero temperature while 
thermodinamically there is no finite-temperature phase transition. 

A. The p-states model 

The simplest way we can introduce interaction among different states in the model is the 
following, consider the quantity 

1 N 

M[{n r }] = -J2(S nr ,o-l/e). (32) 

iV r=l 

Any Hamiltonian of the form 

H = NF(M[{n r }]) (33) 

with F gentle enough, is a good candidate for a mean-field model. We did not try a system- 
atic study of the form ( |33"D for generic F, but we concentrated to the class of monomials, 
where 

H P = -^zi (£(<W - l/e)J (34) 

For p = 1 this model reduces to the BG model. For larger values of p there is interaction 
between different states. The ground state of this model is the same as the one of the 
BG model (all particles occupying the same state) and there are no energy barriers at zero 
temperature. A careful study of the thermodynamics of this model shows that for any p > 1 
there is a first order phase transition from a completely disordered phase with M = for 
T > T c to an 'ordered' phase with M/fl for T < T c . This leads to the curious situation 
that the completely disordered state is dynamically stable at all temperatures but at T = 0. 
This can be understood by a simple argument. Suppose to start the dynamics in a random 
initial condition and consider a sweep that lead to the filling of an empty state. The energy 
variation in this process is 
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(35) 

But according to our hypothesis J2r(^n r ,o — ~) is a random variable of order \fN , correspond- 
ingly 5H ~ N^ip^ 1 )/ 2 and the acceptance rate 

e -m „ e -m^ (36) 

For finite temperature and large N all the moves are accepted and the energy in average never 
decreases. In other words the statistic of configurations is not changed by the dynamics. 
A crossover is found for (3 ~ Nip- 1 )/ 2 showing that the relevant scale of temperature for 
the dynamics is different from that of the statics. Right at zero temperature, where only 
the sign of the energy change and not the magnitude matters, the dynamics of the model 
coincides for any p with the one of the conventional case p = 1. 



B. The effect of metastability 

The p-states model has no metastability at zero temperature. We want to study here a 
simple model where metastability is present but without interaction. In the BG model the 
ground state is reached by emptying progressively more and more states. To empty a given 
state at a certain time t it is necessary to pass in a configuration where a unique particle 
occupies that state. We then consider the following model, 

N 

H = E (-^,o + <?<W) (37) 

r=l 

where g is positive constant and we have the usual constraint eq.(0). At zero temperature 
the transition n r = 2 — > n r = 1 is forbidden, hence energy barriers are present in the model. 
More general models are obtained including in the Hamiltonian all possible terms of the type 

&n r ,k 

N oo N 

H = ~ 9k$n r ,k = ~ 9n r ( 38 ) 

r=l fc=0 r=l 
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We focus here to the case (|37|). The statics of this model is easily solved. We obtain the 
free energy, 

f3f = log(» - log{e z + e p - 1 + z{e~^ - 1)) (39) 

and the fugacity is related to the temperature 4 by the (^-independent relation eq.(^). 
The equilibrium probabilities Pj~ (see eq.(§)) are given by, 

p _ zk ~ X exp(/36 k ,o ~ Pg5 k ,i) , , 

k ~ k\(e* + e-to-l) l4UJ 

The dynamics of this model is expected to be substantially different to that of the BG 
model at least at very low temperatures. Concretely, at zero temperature the ground state 
is the same as for the BG case but there is a large number of metastable states (for instance, 
half states empty and half of the states with two particles). It is easy to show that for 
each value of E between E = —1/2 and the ground state E = —1 there exists a metastable 
configuration with that energy. Then we expect the value of the energy extrapolated to 
infinite time to depend strongly on the initial configuration. In order to minimize the energy 
we have to maximize Po and minimize Pi. While the maximization of Po is a process where 
entropy barriers are dominant (this is the reason why the BG model defined as E = — P 
is interesting) this is not the case for minimizing Pi where entropy barriers are absent. 
Then, independently of the initial configuration we expect that Pi will go to zero in the 
exponentially fast for large times. In these conditions, we do not expect that the adiabatic 
solution of (II) can give a good approximation of the dynamics. This approximation was 
based on the fact that in the BG the surfaces of constant energy are connected, a situation 
which does not hold here. 

However the dynamics of this model can be directly solved as in the BG case. Skipping 
all the details we find for the generating function: 

= (x- 1)(-(1 + A(t))^^ + (1 + 2P 2 (e^ - l))G(x,t) + 

X(t) - 2P 2 (e-^ - l)(e x - P (x - 1)(1 - e-? {1+9) )) (41) 
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where A(i) = Po(^)( e_ ^ 1+9 ^ ~~ !)• Observe that eg. fl4ll) depends only on the probabilities 
Po{t) and Pzit). The solution is more complicated to that of eg . C|X3]) , however it can be 
found. In figure 3 we show the result of the numerical integration of (|4~i~D for the energy 
at T = compared with the Monte Carlo simulations, starting from a completely random 
initial condition. The energy seems to converge to a value lim^oo E(t) = —.564, a result that 
it would be interesting to derive analytically. It is under current study the finite temperature 
dynamics, where we expect the effect of the energy and entropy barriers to combine to give 
rise to a dynamics slower than that of the BG model. 

V. CONCLUSIONS 

In this paper we have derived the exact mean-field eguations of the dynamics of the 
Backgammon model. This has been achieved through the study of the single site occupation 
number probability for which a hierarchical set of eguations hold. With the method of 
the characteristic function, we have derived a closed functional eguation for the energy. 
This, although non-Markovian, has a causal character and can be integrated step by step 
discretizing the time. The non-Markovian character of the evolution eguation, suggests that 
history dependent effects should be observable in the system. However the analysis of II, 
where the evolution of the energy was described by an approximate eguation, shows that 
even in subtle phenomena as hysteresis cycles in cooling-heating processes, history dependent 
effects are very small. This should reflect in the fact that the memory kernels that appear 
in the eguation for the energy are short range in time. 

The method of the generating function also allowed us to derive also a system of linear 
Volterra eguations describing the evolution of the energy autocorrelation function. The 
numerical solution of these eguation confirmed the aging behaviour found in I. It would be 
interesting to derive analytically the scaling Cs(i, i') = fit' /t). 

In the last section we have derived the mean-field theory for a model where entropic and 
energetic barriers are combined. We have seen that at temperature zero, starting from a 
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random configuration the system fails in finding the ground state. For future work it is left 
the study of this model for finite temperature. 

Non linear equations with memory appear in phenomenological glass theory under the 
name of Mode Coupling Theory jl]. Mode coupling equations appear naturally in the mean- 
field treatment of the dynamics of disordered [0,0] or quasi-disorderd systems |T7| , in off- 



equilibrium situation they involve a set of coupled integral equation for the two time auto- 
correlation function and its conjugated response function. The most striking manifestation 
of the importance of memory effects in off-equilibrium mode coupling theory is in the aging 
behaviour of the response function ||i5| , [T6| . 



Structural glasses are generally classified as strong glasses (Arrhenius behavior of the 
relaxation time) or fragile glasses (Vogel-Tamman-Fulcher behavior of the relaxation time). 
In this classification the BG model is a strong glass. Polymer glasses are fragile glasses which 
show strong aging effects in its physical properties 0. It would be desiderable to know from 
experiments if there is a correlation between the fragility of glasses and its aging properties. 
This could shed light on the role of energy barriers in the mechanism of the glass transition. 
We believe that only entropy barriers cannot yield aging effects in the response function. 
In this framework a more detailed study of the BG model with metastability (as presented 
in the last section) at finite temperature could be instructive. In particular, the study of 
the relaxation time as a function of the temperature and the existence of aging due to the 
presence of energy barriers. 
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VI. APPENDIX A 



In this appendix we derive the evolution equation for the probability Pk(t). Define as 
Nk(t) the number of states occupied by k particles (Pk(t) = Nk(t)/N). 
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The processes leading to a variation of N k can be classified in this way: 

• Process A+: arrival of a particle in a state with k — 1 particles. 

• Process A—: departure of a particle from a state with k particles. 

• Process B+: departure of a particle from a state with k + 1 particles. 

• Process B—\ arrival of a particle in a state with k particles. 

Note that in a sweep the processes above are not mutually exclusive, so, for example the con- 
temporary occurency of A+ and B— lead to no variation in N k . At each Monte Carlo sweep 
three independent random variables are extracted: a departure state d with probability 
rid/N, an arrival state a with probability 1/N and an acceptance variable 



In terms of these variables the variation in N k in each process is given by: 

• Process A+: 5„ a , fc -i[l - S kA (l - <5„ d ,i)4, ] 

• Process A-: -S ndjk [l - £„ a ,o4,o + S ki i + S kj iS na>0 S Xi0 \ 

• Process B+: S ndjk+1 [l - S nat0 S Xi0 + S kt0 S nai oS X}0 ] 

• Process B-: S na>k [l - 6 hfl (l - ^,1)^,0] 

The contribution in the different processes can be easily understood, for example in process 
A+ we must have k — 1 particles in the arrival state. If k = 1 and rid > 1 the move imply 
an energy cost, and is accepted only if x = 1. 

Summing all the contribution and averaging over p, a and x we find: 




1 with prob. e " 
with prob. 1 — e~ 



(42) 



(N k (t + St) - N k (t)) = N[P k (t + St) - P k (t)} = — |p = 
(k + l)(P k+1 - P k ) + P fc _! + P (e _/J - - 4,o - £P fc + (k + l)P k+ i) 




(43) 



Very similar considerations lead to ( PD| ) if one restricts the balance equation to a subset 



of the whole space. 
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VII. APPENDIX B 



In this Appendix we obtain the solution of eq.(|Tl|). We perform the change of variables 
(x,t) — > (u,t) where x — 1 — e u+ Jo ds ( 1 + A ( s )) = e u B(t,0). In terms of the new variables 
eq.(p"T|) reads, 

— = e u B(t,0)(G + X) (44) 

where G(u,t) = G(x(u,t),t) This is a linear differential equation which can be readily 
solved 

G(u,t) = e eU Io dsB ^F(u) + e u f ds X(s)B(s, 0)e e " AvB( - v ' 0) (45) 

J o 

where F is an arbitrary function. 

Going back to (x,t) and imposing the initial condition we get eq.(|i~3|). 
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Figure Captions 

Fig. 1 The decay of the energy at zero temperature starting from a completely random con- 
figuration at time t = 0. We compare the numerical solution of (|I5]) (full lines) with 
the Monte Carlo simulations for N = 10 5 and the integration of the adiabatic equation 
of II with the same initial condition. 

Fig. 2 The correlation function at zero temperature as a function of (t — t w )/t w for different 
t w (t w = 10, 30, 100, 300). We take it as a good indication for the t/t w scaling at large 
times. 

Fig. 3 The correlation function at zero temperature for t w = 10 compared with the montecarlo 
data for N = 10 5 . 

Fig. 4 Energy vs. time in the model with energy barriers (theory + Monte Carlo data at 
N = 10 5 ), starting from a random configuration at time zero. We observe exponential 
decay to E = —.564. 
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